Modeling friction on a mesoscale: Master equation for the earthquake-like model 
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The earthquake-like model with a continuous distribution of static thresholds is used to describe 
the properties of solid friction. The evolution of the model is reduced to a master equation which 
can be solved analytically. This approach naturally describes stick-slip and smooth sliding regimes 
of tribological systems within a framework which separates the calculation of the friction force from 
the studies of the properties of the contacts. 
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In spite of its crucial practical importance, friction is 
still not fully understood [l[ . It raises questions at many 
scales, from the atomic scale studied nowadays by atomic 
force microscopy to the macroscopic scale of a solid block 
sliding on an other. A simple mesoscopic model has 
been introduced to bridge the gap in scales and describe 
the main experimental observations, such as stick slip or 
smooth sliding, in terms of the properties of local con- 
tacts. This widely used Burridge-Knopoff spring-block 
model , initially introduced to study earthquakes (EQ 
model), has been developed by Olami, Feder and Chris- 
tensen Q. It describes the contacts in terms of elastic 
springs and junctions that break at a critical force. Com- 
puter simulations 0, [j| showed that the EQ model may 
reproduce the experimentally the observed stick-slip and 
smooth-sliding regimes, including the role of velocity and 
temperature, if the model is at least two-dimensional and 
various assumptions on the properties of the contacts are 
made. 

The drawback of such a simulation approach is that 
heavy calculations with different parameter sets or con- 
tact properties are required to determine the main fea- 
tures of the model, and it is hard to draw conclusions 
of general validity. The calculations may be tedious be- 
cause a large number of contacts and investigations on 
very long evolution times are necessary to get meaning- 
ful statistics and to make sure that the calculation has 
reached asymptotic properties which are not influenced 
by the initial conditions. Moreover almost all the stud- 
ies based on the EQ model assume for simplicity, and to 
reduce the parameter space to explore, that all contacts 
have identical properties. It turns out that, as we show 
below, this limit is singular and may lead to qualitatively 
incorrect conclusions. 

Here we introduce a master equation (ME) approach 
which is much more efficient than simulations and can 
be solved analytically in cases which are particularly rel- 



* obraun@iop.kiev.ua http : //www. iop .kiev.ua/~obraun 
t Michel.Peyrard@ens-lyon.fr 




k 




FIG. 1: (color online): The earthquake-like model of friction. 



evant. It provides a deeper understanding of friction an- 
alyzed at the mesoscale in terms of the statistical prop- 
erties of the contacts. This splits the study of friction in 
two independent parts: (i) the calculation of the friction 
force given by the master equation provided the statisti- 
cal properties of the contacts are known, (ii) the study 
of the contacts and their statistics, which needs inputs 
from the microscopic scale. Many aspects such as the 
interaction between the contacts and their aging can be 
studied separately to determine their role on the statisti- 
cal properties of the contacts and then accounted for by 
the master equation approach. 

Earthquake-like model. The EQ model describes the 
contact interface, i.e. the interface between the bottom 
of the solid block and the fixed substrate (Fig. [T]). It 
assumes that the interaction occurs through N c asperi- 
ties that make contacts with the substrate. Each asper- 
ity is characterized by its contact area Ai and an elastic 
constant fcj , schematized by an elastic spring on Fig. [TJ 
which can be estimated from hi ~ pc 2 ^/A~i, where p is 
the mass density and c is the transverse sound velocity 
of the material which forms the asperity When the 
bottom of the solid block is moved by X, the stretching 
Xi of an asperity, i.e. its elastic deformation with respect 
to its relaxed shape, increases. The force at the contact 
grows as /j = fcjCCj until it reaches the threshold value 
f s i oc Ai at x s i = fsi/ki oc \J~Ai, at this point the contact 
rapidly slides, and /$ and xi drop to a small value before 
a contact is formed again. 

Let Pcixg) be the normalized probability distribution 
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of values of the thresholds x s % at which contacts break. 
The model is studied in the quasi-static limit where iner- 
tia effects are neglected so that only positive values of Xi 
are relevant. The distribution P c (x) can be characterized 
by its average value x s and standard deviation a s . A typi- 
cal example is a Gaussian distribution P c (x) — G(x] x, a) 
0. 

To describe the evolution of the model, we introduce 
the distribution Q(x;X) of the stretchings Xi when the 
bottom of the solid block is at a position X. It is 
normalized by J Q dx Q(x; X) = 1 for all X. Let all 
asperities be initially relaxed or weakly stressed, e.g., 
let the distribution Q(x;0) = <5i n i(cc) be the Gaussian 
Qini(x) = G{x; x iai , a- ini ) with x ini < x a and cr ini < a s . 
Now, let us adiabatically increase the displacement X of 
the bottom of the solid block while the substrate remains 
fixed. The sum of the elastic forces exerted on the bot- 
tom of the block by the stretched asperities makes up the 
friction force 



can be written as 



F(X) = N c (ki) / x Q(x; X)dx . 



(1) 



The evolution of the system, deduced from the numer- 
ical simulation of the EQ model is shown in Fig. [2^. It 
shows that, in the long term, the initial distribution ap- 
proaches a stationary distribution Q s (x) and the total 
force F becomes independent on X. The final distribu- 
tion is independent of the initial one. An elegant mathe- 
matical proof of this statement was presented in Ref. 7] . 
The statement is valid for any distribution P c {x) except 
for the singular case of P c (x) = 5{x — x s ). 

Master equation. Rather than studying the evolution 
of the distribution Q(x\X) by a simulation of the EQ 
model it is possible to describe it analytically. Let us 
consider a small displacement AX of the bottom of the 
solid block. It induces a variation of the stretching Xi of 
the asperities which has the same value AX for all as- 
perities if the deformation of the bottom surface of the 
block can be neglected. As discussed below, the general 
case where the relative positions of the asperities on the 
surface are allowed to vary can be cast into a general- 
ization of this formalism. The displacement X leads to 
three kinds of changes in the distribution Q(x; X): first, 
there is a shift due to the global increase of the stretch- 
ing of the asperities, second, some contacts break because 
the stretching exceeds the maximum that they can stand, 
and third, those broken contacts form again, at a lower 
stretching, after a slip at the scale of the asperities, which 
locally reduces the tension within the corresponding as- 
perities. These three contributions can be written as a 
master equation for Q(x; X): 

Q(x;X+AX) = Q(x-AX; X)-AQ_(x; X)+AQ+(x; X). 

(2) 

The first term in the r.h.s. of Eq. @ is just the shift. 
The second term AQ_(x; X) designates the variation of 
the distribution due to the breaking of some contacts. It 



AQ-(x; X) = P(x)AXQ(x;X) 



(3) 



where P(x) AX is the fraction of contacts that break 
when the position changes from X to X + AX. According 
to the definition of P c (x) the total number of unbroken 
contacts when the stretching of the asperities is equal to 
x is given by N c P c (£.) d£. The contacts that break 
when X increases by AX, so that the stretching of all 
asperities increases by AX, is the number of contacts 
which have their thresholds between x and x + AX , i.e. 
JV c P c (j:)AI. Thus 



p(x)=p c (x)/ r<%p c (o 



(4) 



The broken contacts relax and have to be added to the 
distribution around x ~ 0, leading to the third term in 
Eq. We denote by R(x) the normalized distribution 
of stretchings for the relaxed contacts. Writing that all 
broken contacts described by AQ_(x; X) reappear with 
the distribution we get 



AQ+(x;X) = R(x) d(AQ^;X 



(5) 



Equation (|2|) can be rewritten as [Q(x;X + AX) — 
Q(x; X)} + [Q{x- X) — Q(x — AX; X)} = -AQ_ (x; X) + 
AQ + (x; X). Taking the limit AX — > 0, we finally get the 
integro-differential equation 



dQ{x;X) dQ(x;X) 



dx 



dX 



+ P(x)Q(x;X) 



R(x) / d^P(0Q(tX), (6) 
Jo 



which has to be solved with the initial condition 
Q(x; 0) = Qini(^)- Notice that Qim(x) cannot be an ar- 
bitrary function, because the contacts that exceed their 
stability threshold, must relax from the very beginning. 

Once the distribution Q(x; X) is known, we can calcu- 
late the friction force F(X) from Eq. (P). The static 
friction force is the maximum of F(X), i.e., F s — 
F(X S ), where X s is a solution of the equation F'(X) = 
dF(X)/dX = 0. In order to simplify the calculation, 
we will assume in what follows that R(x) = S(x), i.e. 
when the broken contacts stick again the asperities are 
completely relaxed. 

Analytical solutions of the master equation can be ob- 
tained for some particular cases of contact properties 0] , 
such as a rectangular P c (x) distribution. Moreover, for 
one particular but important choice of the initial dis- 
tribution, when all contacts are relaxed at the begin- 
ning, Qi n i(a;) — S(x), we can find analytically the ini- 
tial part of the solution in a general case [8|. For the 
rather general case of a Gaussian distribution of thresh- 
olds, P c (x) = G(x]x s ,o~ s ), a numerical solution of the 
master equation ^ is presented in Fig. [5b. One can 
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FIG. 2: (a) Evolution of the EQ model. The curves show the distribution Q(x;X) for incrementally increasing values of X 
(with the step AX ~ 1). The distribution P c (x) is Gaussian with x s = 1 and a s = 0.05, the initial distribution Qi n i(aO is 
Gaussian with X[ n i = 0.1 and <jj n i = 0.025. (b) Solution of the master equation with the increment AX = 1.09 for the same 
model parameters. The bottom panels show the corresponding dependences F(X) for {ki) = 1 and K — do. 
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FIG. 3: (color online): The final distribution Q(x) for the 
parameters from Fig. [2] (solid curve; crosses show the averaged 
final distribution for the EQ model). The red dotted curve 
shows the distribution P c {x), and the blue broken curve shows 
P(x). 

see that it is almost identical to that of the EQ model 
(Fig. [2^i), except for the noise on the EQ distributions. 
The distribution Q{x\X) always approaches a station- 
ary distribution Q s (x). The final distributions of the EQ 
model and the master equation approach are compared 
in Fig. rj 

The steady-state, or smooth-sliding solution, i.e. the 
solution of Eq. j[6]) which does not depend on X, can 
easily be found [J] . In a general case it can be expressed 
as 

Qs(x) = C- 1 e(x)E P (x), (7) 

where Q(x) is the Heaviside step function, Ep(x) = 
e~ u i x ), U(x) = f*d£P(Z), and C = J™dxE P (x). The 



(kinetic) friction force is then equal to 

/>oo 

F={N C /C){h) / dxxE P {x). (8) 
Jo 

In the general case, let the distribution P c {x) be of 
bell- like shape with the maximum at x s and the width a s . 
When X shifts for the distance x s , due to the breaking 
and reforming of contacts with a lower stretching, an ini- 
tially peaked distribution Q(x; X) broadens by the value 
~ <j s (Fig. [5]). Therefore, any initial distribution tends to 
the stationary one as \Q(x;X) — Q s (x)\ oc exp(— X/X*), 
where X* ~ x 2 s /a s . 

Thus, in a general case the solution of the master 
equation always approaches the smooth-sliding one given 
by Eq. ([7]). However, there is one exception from this 
general scenario. When all contacts are identical, i.e., 
all contacts are characterized by the same threshold x s , 
P c (x) = S(x — x s ), the model admits a periodic solution 
8]. This singular periodic solution has been found in 
simulations and analyzed as describing the stick slip [9(, 
but actually it is unphysical and ceases to exist as soon 
as non-equivalent contacts are considered, whatever their 
precise properties. As discussed below the stick slip can 
be deduced from the solution of the master equation, but 
its origin is different [lfj. 

Stick slip and smooth sliding. The master equation 
allows us to compute the friction force F(X) when the 
bottom of the solid block is displaced by X. But actually 
we don't control X. The displacement is caused by a 
shearing force Ft applied on the top of the solid block 
which displaces the top surface by Xt- As the strain 
on the solid is generally small the deformation of the 
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solid block can be assumed to be elastic so that Xt is 
related to the applied force by Ft = K(Xt — X) where 
K is the shear elastic constant of the solid block. The 
total force applied to the bottom of the solid block, which 
determines its displacement X, is the sum of the applied 
force and the friction force 

F tot = K(X T -X)- F{X) . (9) 

It can be viewed as deriving from the potential 

V{X T ,X)= l -K{X T -X) 2 + J F(£)dt, (10) 

which determines the behavior of the solid block sub- 
jected to friction and applied force. A necessary condi- 
tion for smooth sliding is that Xt and X grow together 
with Xt — X = B, where B is a constant that mea- 
sures the shear strain of the solid block during the slid- 
ing. It is determined by the condition dV /8(Xt — X) = 
—dV/dX = , which simply means that the total force 
on the interface vanishes. Smooth sliding also requires 
this state be stable, 

d 2 V d 2 V 

d[ x T -xy = dx^ Q OT **W>-K- (ID 

If we start from relaxed asperities, in the early stage of 
the motion F(X) is a growing function of X, and then it 
passes by a maximum when some contacts start to break 
and reform at lower asperity stress. As a result F'(X) 
becomes negative. Depending on the value oi K two situ- 
ations are possible. For large K (stiff block) F'(X) never 
falls below —K and the smooth sliding is a stable steady 
state. For small K (soft block) F'(X) can become smaller 
than —K so that the stability condition [[lip is no longer 
valid. The instability causes Xt ~ X to change abruptly 
by a breaking of all the contacts and a quick slip of the 
block before the contacts reform with relaxed asperities. 
And the process can repeat again, leading to the famil- 
iar stick slip motion. The master equation, which gives 
F{X) can be used to compute the period of the stick 
slip, and, when the asperities fully relax before the con- 
tacts reform, an analytical solution can be obtained Q. 
It should be noticed that the existence of a stick slip is 
not only determined by the stiffness K of the solid block. 
The properties of the asperities, defined by the distribu- 
tion P c (x) of the stretching for which they break is also 
essential because it determines the expression of F(X) 
and hence the minimal value of F'(X). 

Discussion. The ME formalism can be extended 
to take into account various generalizations of the EQ 
model. For instance, in establishing the master equation 
we assumed that the asperities were fixed with respect to 
the solid block and were moving together with it. This 
is an approximation and one can, in principle, take into 



account the elastic deformation of the interface by in- 
troducing for instance a position dependent distribution 
Q[x,X(r)] where r denotes the position of an asperity 
on the interface, and X(r) the local translation of the 
interface averaged over a mesoscopic scale. The master 
equation must then be coupled to an equation describ- 
ing the elastic deformation of the interface, subjected to 
the contact forces at each point r. This illustrates the 
new viewpoint introduced by the ME approach which de- 
scribes the phenomena at an intermediate scale between 
the microscopic scales of the contacts and the macro- 
scopic scale of the displacement of the solid block. The 
ME equation only has a meaning if, on this intermediate 
scale, there are many individual contacts, allowing us to 
study them as a statistical distribution and not individu- 
ally. A simpler view of the effect of the elastic interaction 
between the contacts is to describe it as a renormalization 
of the distribution P c (x) in a mean field like approach. 

Another aspect which can be introduced in the ME 
formalism is the aging of the contacts [8| . Experiments 
and MD simulations show that the static friction force 
grows with time since a contact is formed. As a result, 
if newly formed contacts are characterized by a distribu- 
tion P c i(x) this distribution evolves with time, moving 
to a higher mean value and a smaller standard devia- 
tion. Aging is a stochastic process which can be described 
by a Smoluchowsky equation. The master equation for 
Q(x,X) must then be completed by an equation for the 
evolution of P c {x), which in turns affects P(x) in the 
master equation. Time enters through the velocity of 
the sliding, a faster sliding giving less time for the aging 
of the contacts. The full development is too long to be 
given here [8J] but it is easy to realize that a larger veloc- 
ity leads to a smaller friction force because the contacts 
have less time to age, which is the source of a poten- 
tial instability. Aging is therefore another cause for the 
stick-slip behavior. 

The ME formalism can also accommodate tempera- 
ture effects which enter again through their effect on the 
distribution P(x) because thermal fluctuations allow an 
activated breaking of contacts for asperities which are 
still below their thresholds Q. The main point that we 
would like to stress is that this formalism introduces a 
new viewpoint on the mescoscopic modeling of friction. 
Moreover it splits the analysis of friction phenomena into 
problems that can be studied separately, the statistical 
properties of the contacts, and the evolution of the dis- 
tribution Q which is described by the master equation, 
coupled to additional equations representing different ef- 
fects such as the elastic interactions between asperities, 
the aging of the contacts or temperature fluctuations. 
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